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ion collisions and find them to be quite sensitive to the effects of the medium and the 
production mechanisms. In order to put this on a quantitative footing, as a first step, 
we analyze the azimuthal, transverse momentum, and rapidity correlations of heavy 
quark- anti quark {QQ) pairs in pp collisions at 0{al). This sets the stage for the 
identification and study of medium modification of similar correlations in relativistic 
collision of heavy nuclei at the Large Hadron Collider. Next we study the additional 
production of charm quarks in heavy ion collisions due to multiple scatterings, viz., 
jet-jet collisions, jet-thermal collisions, and thermal interactions. We find that these 
give rise to azimuthal correlations which are quite different from those arising from 
prompt initial production at leading order and at next to leading order. 
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1. Introduction 

The study of relativistic heavy ion collisions and quark gluon plasma (QGP) is 
approaching its zenith with the first experiments performed at the Large Hadron 
Collider at CERN Geneva (though not yet at the top energy) involving lead nuclei. 
Together with the wealth of data already accumulated at the Relativistic Heavy Ion 
Collider at Brookhaven National Laboratory, we now have an enormous task to decipher, 
analyze, and quantitatively explain these observations and extract information about 
the properties of the QGP. These analyses are also paving the way for additional 
measurements, some of which can already be performed using the present detector 
set-ups, while others will become amenable to studies with the upgrades planned for all 
the major experiments, ALICE, PHENIX, and STAR, etc. Taken in its entirety, this 
represents the most important and fruitful international collaboration in high energy 
nuclear physics to date. 

The focus has now progressed from models to theories and from qualitative to 
quantitative determination of various properties of quark gluon plasma. Enormous 
strides made towards exploring the shear viscosity [1] of the matter produced in these 
collisions is one such example. 

In the present work we consider using heavy quarks to probe the QGP. The 
heavy quarks (only charm and bottom quarks are considered here) offer several unique 
advantages. The conservation of flavour in strong interaction dictates that they are 
produced in pairs {QQ)- Their large mass provides that necessary for their 
production is large and thus one may confidently use pQCD, for these studies. Their 
large masses ensure that the hadrons containing the heavy quarks stand out in the 
swarm of pions. 

Their large mass also provides that, even though buffeted by light quarks and gluons 
during their passage through the quark gluon plasma, the direction of their motion may 
not change substantially. This should make them a valuable probe for the properties 
of the plasma which depend on the reaction plane. Our understanding of the effect 
of the dead cone [21 [Sj H] on the suppression of radiation has undergone quite some 
evolution since it was proposed earlier. It is also not yet clearly established that heavy 
quarks will completely thermalize in the plasma formed at RHIC and LHC energies 
(see Ref. [5]). However it must be safe to assume that the drag [6] suffered by heavy 
quarks will mostly slow it down and the so-called diffusion [7] processes will not alter 
its direction considerably. Thus, the azimuthal correlation of heavy quarks integrated 
over pt may be reasonably immune to the energy loss suffered by them. 

The heavy quarks could be influenced by the flow [8] generated in such collisions. 
If this is true, then a very interesting situation may arise for heavy quarks which is not 
possible for light quarks or gluons. Consider a QQ pair produced in a central collision 
having y = 0. At leading order, their transverse momenta would be equal in magnitude 
and point towards opposite directions. Consider a heavy quark Q moving away from 
the centre with momentum px. Then its partner Q would move with momentum — Pt 
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Figure 1. Energy dependence of the charm quark production in pp collisions. 



towards the centre. Their velocities would be vq = ±Pt /Mt, where Mt = yPx + Mq, 
and Mq is the mass of the heavy quark. Let the radial flow velocity be Vf. Now if 
|vf I > |vq|, the Q will turn back and start moving away from the centre! Thus the QQ 
pair, which should have appeared back-to-back would appear as moving in the same 
direction. This would drastically alter the azimuthal correlation of the pair. A similar 
change of direction of motion is not possible for light quarks and gluons as they move 
with the speed of light. Taking, for example, |vf | ^ 0.6, (see Ref. [8]) one can see that the 
azimuthal correlation of charm quarks for < 1.2 GeV and for bottom quarks having 
Pt < 3.5 GeV could be considerably modifled from their primordial value. Recalling 
that non-back-to-back heavy quarks are produced from NLO processes (see later also), 
this would introduce an interesting richness in these studies. 

Now consider charm quarks (say) produced from the primary processes gg — )■ QQ at 
leading order and gg — gQQ at next-to-leading order. In the absence of any intrinsic kr 
for partons, the quarks from the flrst process will be produced back-to-back, while those 
from the second process will be mostly collinear and will additionally be accompanied 
with a recoiling parton. A comparison of the energy loss suffered by the recoiling 
parton and the heavy-quarks will allow us to obtain flavour dependence of the energy 
loss. A considerable richness to this picture is added by the realization that the splitting 
g —7- QQ, would produce collinear heavy quarks, while the process gg — )■ QQg, where a 
gluon is radiated by one of the heavy quarks will essentially give rise to a flat azimuthal 
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Figure 2. Transverse momentum distribution of the heavy quarks in the central 
rapidity region in pp collisions at yG~2.76 TeV, for different structure functions using 
NLO pQCD. 




Pt (GeV) p^. (GeV) 

Figure 3. (left) Transverse momentum distribution of D'^-mesons and (right) of D+- 
mesons, in pp collisions for y/s= 2.76 TeV. 

correlation. 

So far we have discussed only the azimuthal correlation of the heavy quarks. A study 
of the transverse momentum of the pair and the rapidity-difference of the pair can help 
us disentangle the LO and the NLO processes. Recall that the transverse momentum of 
the QQ pair would be identically zero at LO and equal to that of the recoiling parton at 
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Figure 4. Transverse momentum distribution of single electrons from pp collisions at 
y/s = 7 TeV. 
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Figure 5. Transverse momentum (left panel) and rapidity distribution (right panel) 
of J/ V' from pp collision at y/s= 7 TeV, using color evaporation model. 



NLO. Deviations from the results for pp collisions at the corresponding centre of mass 
energy in nuclear collisions will provide a measure of medium modifications as usual. 

The entire discussion so far assumes that there may be no additional production of 
heavy quarks after the initial prompt production. Often this production for pp collisions 
is taken as a baseline for the study of nuclear modification factor Raa- It is obvious that 
any additional production of heavy quarks, for example due to multiple scattering of 
high momentum quarks and gluons produced similarly, see Refs. [9l [TOl [TTl [12] or due to 
passage of a high energy quarks or gluons through the QGP [TOlIlS], or due to scattering 
among the thermalized partons, if the temperature is sufficiently large flU\ [TT| [T3] , will 
necessitate revision of our estimates for the energy loss suffered by heavy quarks as they 
traverse the QGP, obtained by analyzing the nuclear modification function Raa, (see 
Ref. P). 

By now there is also a growing realization that Raa is not able to seriously 
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Ai^ (radians) A(|) (radians) 



Figure 6. Azimuthal correlation of charm (left panel) and bottom (right panel) 
quarks at 2.76, 5.5 and 7 TeV for pp collisions. The symbols give the LO values 
for da/d(p — (T-^Q/S{A(j)) where 6{A(p) is the size of bin. 

discriminate between different meclianisms of energy loss and evolution of the system [T5] 
and the correlation of the leading hadrons are slowly emerging as more discerning 
probes [16]. Consider a simple example. We need to know the transverse momentum 
of heavy quarks in pp collisions in order to have a base-line to estimate the nuclear 
modifications. The NLO pQCD results for these are easily approximated hj a. K factor 
multiplying the results for LO pQCD (see eg. Ref. [IT]). Now consider the azimuthal 
correlations of heavy quarks produced in similar collisions. As we discussed above, the 
LO pQCD results for the correlation is a delta function around A(f) = tt. However, we 
shall see that the correlation function estimated at NLO, though still peaking at A(f) = n 
fills up the phase-space from zero to vr with an interesting catenary like structure. 
Considering that one uses deviations from pp collisions to obtain results for nuclear 
modifications, these will have to be quantitatively understood for pp collisions before 
we can accurately decipher the later. 

The present work aims at investigating azimuthal, momentum, and rapidity 
correlations for heavy quark-anti quark pairs for pp collisions and setting the stage 
for the study of the deviations in these due to medium modifications in heavy ion 
collisions at the corresponding energies. We also discuss the complexities arising from 
the additional production of heavy quarks due to multiple scatterings. 

The paper is organized as follows. In the next section we discuss various correlations 
for pp collisions using NLO pQCD. In Sect. 3 we discuss the azimuthal correlations in 
Pb+Pb collisions due to initial production and various multiple collisions. Our results 
for pp and Pb-|-Pb collisions are discussed in Sect. 4 followed by conclusion in Sect. 5. 
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2. Proton Proton Collisions 

The results for particle and photon productions in pp collisions serve as a baseline in 
search for quark-gluon-plasma and other medium effects at the corresponding centre 
of mass energy /nucleon for collision of heavy nuclei. This paradigm may have to 
be modified if the recent suggestions for formation of QGP (perhaps only in high 
multiplicity events), Ref. [18] in pp collisions turn out to be valid. An alternative 
criterion of comparing results for peripheral collisions to those for central collisions has 
also been used with considerable success, with the understanding that the peripheral 
collisions may be considered as a superposition of pp collisions. 

The correlation of heavy quarks produced in pp collisions is defined in general as: 

^^^"^ d?pi(Pp2 dyidy2(PpTid?'PT2 ^ ' 
where yi and y2 are the rapidities of heavy quark and anti-quark and pTi are their 
transverse momenta. 

At the leading order, the differential cross-section for the charm correlation from 
proton-proton collision can be written as: 

= d^pXdy/ ^''''' + P^^^ ■ 



In the above Pxi = Pt2 = Pt and 



dyidy2dpT 



dt 



daij{s,i,u) 
di 

/(I + 5..), (3) 



where Xa and Xh are the fractions of the momenta carried by the partons from their 

interacting parent hadrons. These are given by 

Mr Mr 
Xa = ^(e^i +ey^)- x, = -^{e-y^ + e"^^) . (4) 
y s y s 



where Mt is the transverse mass, ^Jmq+p^, of the produced heavy quark. The 
subscripts i and j denote the interacting partons, and fi and fj are the partonic 
distribution functions for the nucleons. We shall use CTEQ5M structure function, 
though we have checked that similar results are obtained for other modern structure 
functions (see later). The differential cross-section for partonic interactions, daij/dt is 
given by 

da^j _ \Mf 



dt IGtcs"^ 



(5) 



where \Mf is the invariant amplitude for different sub-processes as obtained from 
Ref. [19]. The physical sub-processes included for the leading order, O (a^) production 
of heavy quarks are: 

g + g^Q + Q 

q + q->Q + Q . (6) 
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At next-to-leading order, O (af) subprocesses included are as follows 

g + g-^Q + Q + g 
q + q-^Q + Q + g 
g + q{q) ^ Q + Q + q{q) ■ 



(7) 



We show our results for azimuthal correlation C(A0), where A(/)=|0i — as 
well as rapidity correlations, C{Ay), where Ay=yi — y2, of produced heavy quarks. 
We also present {Arj, A0) correlations in the jet radius parameter, R, where R = 



y {AriY + {A(py along with the transverse momentum, invariant mass, and rapidity 
of the pair. 

We verify the accuracy of our results by evaluating the production of J/ip and 
charm measured recently. 

3. Lead Lead Collisions 

Let us now move towards Pb+Pb collisions under study at the Large Hadron Collider 
(LHC). We have discussed that most of the heavy-quarks and so also quarks and gluons 
having large transverse momenta are produced in initial hard collisions. At the energies 
reached at the LHC, the sheer number of quarks and gluons produced in these collisions 
leads to vehement multiple collisions and gluon multiplication. This then leads to a 
quark-gluon plasma at a very large initial temperature. 

As discussed earlier, we would like to know if these initial temperatures are large 
enough to produce heavy quarks as well (see eg. Ref. [13]). The multiple collisions 
among the very high momentum quarks and gluons (the so called jet-jet collisions) have 
been seen earlier to produce substantial number of heavy quarks. These jets, produced 
at very early times r ^ l/px will have to necessarily pass through the QGP which will 
be formed only after r ^ 0.1 fm/c. Do these lead to a substantial production of heavy 
quarks? Some of these questions have been addressed earlier [9l [TOl [TTl [12] . 

Since those early studies, several new developments have taken place. We now 
know the particle rapidly density (see eg. Ref. [2D]), important to calculate the initial 
conditions, for which only values estimated by several authors were known earlier. There 
has, now, been a growing realization that jet-quenching measured in terms of the nuclear 
modification function Rp^p^ is not able to seriously discriminate between various theories 
of evolution of the plasma and the mechanism of energy loss. Thus correlations are being 
studied more closely to help us in this enterprise. 

Thus we extend our earlier study [lOj to explore the correlations of heavy quarks 
in collision of heavy nuclei due to initial production and various multiple collisions, e.g., 
jet-jet interactions, jet-plasma interactions, and the scattering of thermal partons, to 
see if these processes make large contribution to the correlation. This is important as, 
at least the jet-jet collision was found to make a large contribution to the production of 
heavy quarks [TO] . 
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3.1. Prompt Interactions 

The basic formulation which gives the correlation of produced heavy quarks from initial 
fusion of gluons and quark-anti quark annihilation in proton-proton collision is given by 
Eq. [H Thus the azimuthal distribution of heavy quark for Pb+Pb collision at 6 = is 
given by 

^ d^Pl Cpp2 ^^12 ^3^^ ^3^^ ( ) 

For central collisions of lead nuclei, the nuclear thickness function is taken as Taa= 
292 fm~^. In the above pi and p2 are the momenta of the heavy quarks produced. 

3.2. Jet- Jet Interaction 

The initial hard scattering will produce massless gluons and light quarks in large 
numbers. These partons have large transverse momenta. These quarks and gluons 
may ultimately thermalize because of frequent interactions among themselves and if 
sufficient energy is available, their interactions may lead to the production of heavy 
quarks as well. Here we give the formulation for azimuthal distribution of produced 
heavy quarks pair from jet-jet interaction. Since the jet-jet contribution to the heavy 
quark production is comparable to that of primary production [HI HO], it should be 
interesting to see if their azimuthal distributions differ. 

As a first step we obtain the distribution of light partons, having px > 2 GeV, from 
a LO pQCD calculation using CTEQ5M structure function, for pp collisions at 2.76 TeV 
and 5.5 TeV. We parametrize them as: 



dN 



dyd^PT d'^prdy 



y=0 



[i + ptIby 

= h-Avr) , (9) 

where the K factor is taken as 2.5 to account for higher order effects and the parameters 
C, 5, and /3 are given in Table 1. The factorization and renormalization scales are chosen 
as Q = pt 

Now the azimuthal distribution of heavy quarks for collisions having an impact 
parameter, 6 = 0, due to jet-jet interaction can be written as: 

dN 1 /■ „,4_„ f d^Pad^Pb^A, 



a^Pid'^P2 16(27r)* J J uja^b 



X 



9lflt{PTa)flApi 



+ 9lT.[flk^PTa)fl;,^{PTb) M,,^QQ V (g, ^ g,) 



2 



(10) 
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Figure 7. Transverse momentum, invariant mass and rapidity distribution of charm 
and bottom quark pairs at LHC. 



where Pa, Pb are the four momenta of the incoming partons and pi and p2 are the same 
for the outgoing heavy quarks, and stands for the flavour of the hght quarks. The jet 
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Figure 8. Rapidity difference of Ay = Uq — Uq, of charm and bottom quarks at LO 
and NLO in pp collisions. 



distribution function fj^^ipx) is given by 

(27r)3 dNi 



fUPi 



'Ill 



This follows the Bjorken space-time correlation used earlier in Refs. |9l [TOl |2T]. In the 
above, pt is the transverse momentum, y is rapidity, t] is the space-time rapidity, and 
Qi is the spin-colour degeneracy of the partons, which is 2 x 8 for gluons and 2x3 for 
quarks. As indicated earlier, Tj and Tf are the the formation time of the jet and the life 
time of the plasma. R'j' is the transverse size of the system. Now the Eq. [TD] reduces to: 



El E2 



dN 



d^Pid^P2 16(27r 
1 



d^x J d^pTh dy\ 



X 



^9gf^et^PTa)f^^^{PTb) 

E{ijtt(^-)^tt(^-) 



M, 



gg->-QQ 



M. 



qq^QQ 



+ 9aX 



+ ill ^ Qi) 



(12) 
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Figure 9. (A77, A0) correlations of heavy quarks produced in pp collisions at ^/s=2.76 
TeV and 5.5 TeV at LO and NLO. 



where d'^x=TdT rdr dr] d(f)r and d^p=pTdpT dcpp Edy. Further 

6{j:e) 



KPTb-PTbo) 



[Pti cos(0i - + pt2 cos(02 " 4>b) " Mn cosh(?/i - r]) - Mt2 cosh(?/2 - v)] 

(13) 

Thus the final integration obtained by putting the above expression in Eq. 
reduces to: 



El E2 



dN 



ln(r//i 



d^Pid^P2 16(27r)27ri?2 



dr] ( 



X 



PTa[pTi cos(0i - (pb) + Pt2 cos((/)2 - (j)b) - Mti cosh(?/i - rj) - Mt2 cosh(j/2 - ri)] 



2^|et(P^«)^|et(P™) 



M. 



gg^QQ 
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Figure 10. (Colour on-line) Azimuthal correlation of heavy quarks from prompt 
interaction for lead on lead collisions at LHC, having different transverse momenta 
and rapidities close to zero. The symbols give the corresponding LO values, with the 
same bin-size for A0. The upper panels are for 2.76 ATeV while the lower panels 
give results for 5.5 ATeV. The left panels give results for charm quarks while the right 
panels give the results for bottom quarks. 



+ {qi ^ Qi 



with pTbo given by, 

VTbQ = 



(14) 



Pti cos( 



PtiPt2 cos(0i - 02) - MtiMt2 cosh(?/i - y2) - m^Q 

- 4>b) + Pt2 cos(02 - 06) - Mti cosh(?/i - rj) - Mt2 cosh(?/2 - rj) 

(15) 



These are similar to the expressions obtained earher by Levai et al [TT]. 

The formation time for the jets (hght pt partons) is taken as Tj=0.1 fm/c, as we 
count those having > 2 GeV, as jets. We take tj ^ _Rrj, of the system and perform 
rest of the integration numerically. Note, however, that Tj and Tf appear in the above 
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Figure 11. Aziniuthal correlation of Ireavy quarks from jet-jet interaction for lead on 
lead collisions at LHC, for different transverse momenta. 



expression only through the term ln(rj/rj) and thus reasonable variations in their values 
will lead to only a modest variation in the results. 

3.3. Jet-Thermal Interaction 

Now, we consider passage of high energy energy jets through quark gluon plasma and 
estimate azimuthal dependence of the produced heavy quarks. 

Suppose, a light parton is produced at position r, and moves at an angle a where 
cos a=f.d, then the distance d travelled by the jet along the direction d, before it reaches 
the surface is given by: 

d = —r cos a + \J B? + r'^ sin^ a , (16) 

Hence the time available for heavy quark production is the minimum of [r^^, rj], where 
Td is time taken by the parton to cover the distance d and tj is the time when quark 
gluon plasma hadronizes. 
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Figure 12. Azimuthal correlation of Ireavy quarks from jet-tliermal interaction for 
lead on lead collisions at LHC, for different transverse momenta. 



The azimuthal distribution of the produced heavy quark from jet-thermal 
interaction is given by 

dN 1 



El E2 



d?pi(Pp2 16(27r) 



d X 



hdyb 

PTbO 



[Pti cos(0i - 0b) + Pt2 cos(02 - 0b) - Mti cosh(?/i - yb) - Mt2 cosh(?/2 - Vb)] 



299f^et^PTa)f^l,{PTto) 
+ 9lT.{f^et^PTa)r;i{PTm 



M. 



gg^QQ 



M, 



qq-^QQ 



(17) 



The pt distribution for jet partons is given in Eq. O The distribution of the thermal 
partons and the cooling of the plasma is given in the Sect. 3.4. 
After some simplifications, the final result is given by 
dN 1 



El E2 



d^Pid^P2 1Q{2tt)HR^ 



rdr dr drj d(f)f, dy^ 



Correlations of. 



16 





Figure 13. Azimuthal correlation of Ireavy quarks from ttiermal interaction for lead 
on lead collisions at LHC, for different transverse momenta. 



X 



PTa[PTl COS(01 



PTW 

b) + Pt2 cos((/)2 - (pb) - Mti cosh(?/i - yh) - Mt2 cosh(?/2 - Vb)] 

2 



99- 



+9,j:[h%t(PTa)f!liPTb0) 

which is then evaluated numerically. 



:i8) 



3.4- Thermal Interaction 

We have discussed earlier that the multiple scatterings among the quarks and gluons 
leads to the formation of quark gluon plasma at a large initial temperature. An 
interaction among the thermalized partons may also lead to charm production provided 
the initial temperature of quark gluon plasma is high. Using the recent results from 
ALICE at ^/s=2.^Q A TeV for central collisions of lead- lead nuclei, we take particle 
multiplicity density to be dN/dy=2850 at ^/s=2.76 TeV/nucleon, [20] and extrapolate 
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Figure 14. Relative contributions of prompt, jet-jet, jet-thermal, and thermal 
productions of charm quarks having transverse momenta of about 2 GeV/c (left panel) 
and 4 GeV/c (right panel) to their azimuthal correlation. 



it to 3000 for ^/s = 5.5 TeV/nucleon. Now using the relation [22] 

— = 4aToVo (19) 



45C(3)7r/2^ dy 

and initial formation time for QGP, ri=0.1 fm/c, we estimate To to be 653 MeV at 2.76 
TeV/nucleon and 664 MeV at 5.5 TeV/nucleon respectively. 

Recall also that at RHIC energies, Xj up to 0.6 fm/c have been used, specially to 
for the part of the evolution which could be described using hydrodynamics. One may 
imagine Xj getting smaller at LHC energies, due to increased activity of minijets, etc. 
Thus for example, the parton saturation models [23] suggest that pgg^^ at LHC energies is 
close to 2 GeV, which suggests that the initial time r, for the plasma would be ~ 1/Psat 



PbPb@LHC, Vs=2.76 A TeV 
lyl < 0.5, < 4.0 GeV 
p„=Y|3,ni^ 
y=(l-|' 

= P^x(r/Rj)' 
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Figure 15. Azimuthal correlation of heavy quarks from prompt interaction for lead 
on lead collisions at LHC, with flow 
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Table 1. Parametrization of the quark and gluon distributions from initial scattering 
of partons at 2.76 and 5.5 TcV in pp collisions, using LO pQCD and CTEQ5M structure 
functons, for px > 2 GeV. 







C [l/GeV2] 


B [GeV] 


/3 




u 


1.078x103 


1.127 


5.615 




d 


1.279x10^ 


1.099 


5.579 


2.76 


s 


1.395x10^ 


1.899 


6.432 




u 


3.371x102 


1.434 


5.999 




d 


3.734x10^ 


1.401 


5.953 




o 

to 


2.917x10'^ 


1.892 


6.523 




u 


7.961 xl02 


1.293 


5.580 




d 


9.478x10^ 


1.254 


5.539 


5.5 


s 


1.228x10^ 


2.174 


6.418 




u 


2.659x102 


1.663 


5.966 




d 


2.908x10^ 


1.624 


5.924 




g 


2.449x10^ 


2.192 


6.519 



or about 0.1 fm/c. We shall discuss the consequences of taking large formation times 
(see later). 

Thus the azimuthal distribution of heavy quarks produced from interactions of 
thermalized partons is given by 

dN 1 



El E-2 



d^Pid^P2 16(27r)8 



d X 



x 



hdvb 

PTbO 



\Pti cos(0i - (/)b) + pt2 cos(02 - 4>b) - Mti cosh(|/i - Ub) - Mt2 cosh(|/2 - Vb)] 



X 



M. 



gg^QQ 



M, 



qq- 



{<li ^ Qi) 



+ 9'qT.{ftliPTa)fl{PTb0) 
i 

where (the boosted) thermal distribution of partons is approximated as 

fthiPT^y^v) = exp 



(20) 



~ cosh(y - r]) 



(21) 



The above integration is done numerically, with the temperature varying according 
to Bjorken's cooling law, i.e. T^r = constant, till the temperature drops to about 160 
MeV. 
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4. Results 



4.1. Proton Proton Collisions 

In the results to be reported in the following, we shall use the CTEQ5M structure 
function, though some results are also given for other structure functions. The mass of 
the charm quarks is kept fixed at rric =1.5 GeV, while that for bottom quarks is mi, = 
4.5 GeV. The factorization and renormalization scales are taken as C ^ rrtq + p'^ with 
factor C = 2 for charm quarks and 1 for bottom quarks. The NLO pQCD code (NLO- 
MNR) developed by Mangano et al. [211 125] has been used for the initial production of 
heavy quarks. 



4-1.1. Production of heavy quarks, charmed mesons, and J /if): The results for charm 
production along with recent results obtained at LHC for pp collisions are shown in 
Fig. [H For the sake of exploration we have also included results for rUc =1.2 GeV and 
the structure function CTEQ5M. A very good description of the data Ref. [26], without 
any adjustment of parameters is seen (see also Refs. [TTj l27]). 

We have given the results of our calculations using several structure functions in 
Fig.[2]for the production of charm and bottom quarks at central rapidities in pp collisions 
at 2.76 TeV. We see that use of any of the more modern structure functions gives results 
which differ by just a few percent from each other. 

One may also consider the production of D-mesons by writing schematically: 

£|^ = B«^®fl(Q^J/«). (22) 

where the fragmentation of the heavy quark Q into the heavy-meson Hq is described by 
the function D. We have assumed that the shape of D{z), where z = Pd/Pc, is identical 
for all the D-mesons 



z[l-l/z-ep/{l- z)]^ 



€„ is the Peterson parameter and 



dzD{z) = l. (24) 

The production of a particular D-meson is then obtained by using the fraction for it, 
determined experimentally [291 EO] . 

A comparison of our results for Z)" and production with the preliminary data 
obtained by ALICE experiment [31] is shown in Fig. |3l We give results for = 0.001, 
0.06, and 0.12 to show the sensitivity of our calculations to this variation. Considering 
that no parameters have been adjusted, the results seem to be satisfactory. More detailed 
and accurate data will definitely put stringent constraints on all the inputs. 

Note that the semi-leptonic decay of D-mesons has been extensively used to study 
the production of charm and bottom quarks, as well as the energy loss suffered by them. 
The electrons coming from charm decay, for example, are obtained by convoluting the 
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distribution of D-mesons (Eq. [22]) with the electron decay spectrum [22] and accounting 
for the branching to a particular D-meson [291 130] . In case the contributions of the B and 
the D mesons can not be distinguished, one should use the B and D-meson mixtures, 
with appropriate branchings, B ^ e, D ^ e and B ^ D ^ e. The semileptonic 
decay of B-mesons becomes important at higher pt in spite of their reduced production, 
though the contribution of the B ^ D ^ e channel drops rapidly with increase in pT 
(see e.g., Ref. [33]). 

The ALICE experiment has, however, obtained the single electrons from the process 
c — 7- — e [34j. The upgrades of STAR and PHENIX experiments at RHIC will also 
be able to measure this. 

In Fig. m we compare our results for the electrons measured by the ALICE 
experiment with the decay of charm and a reasonable agreement is seen. In a future 
publication, we shall report on the consequences of introducing an intrinsic kr for the 
partons and also using different parametrizations of the decay spectrum of the electrons. 

The production of J/ in pp collisions is yet another important observable, which 
is closely related to the production of charm quarks. For example, using the colour 
evaporation model, one can write: 



where M is the invariant mass of the pair, y is its rapidity, mn is the mass of D- 
mesons, and F is the (constant) colour-evaporation factor which should be fixed by 
evaluation at some energy. There is one small detail which should be mentioned here; 
the LO pQCD calculations for heavy quark production produce cc pairs with pair- 
momentum identically equal to zero (though the NLO processes do provide them with a 
net transverse momentum). This is corrected by imparting an intrinsic fc^ to the partons 
(see e.g. [35]). Only for these calculations we impart an intrinsic /ct of 1.5 GeV/c to the 
partons. 

We show our results for the transverse momentum and the rapidity distribution of 
J/ip in Fig. [5] along with the experimental results for pp collision at 7 TeV obtained 
for prompt J/tp by the LHCb experiment [36]. (Note that the ALICE collaboration 
has measured the inclusive J/ip which includes the b-decays [37]. Even though this 
contribution is of the order of 10%, it is often accounted for by adding the b — )■ J/ip 
contribution measured by the LHCb experiment.) We have explored the consequences of 
varying the intrinsic on the pt distribution of J/V' and as expected the slope of the pr 
distribution decreases with increase in fc^. A reasonable description of the distribution of 
the transverse momentum and the rapidity distribution is seen. An accurate description 
of the data will involve a more detailed exploration of the parameters. For example, the 
colour evaporation coefficient is kept fixed in these calculations, to magnify the effect 
of varying intrinsic /ct- Of-course the change of intrinsic /ct will not affect the rapidity 
distribution. 

It will be interest to continue with this study for the prompt production of higher 
resonances of cc as well as of bb, when more accurate and detailed data become available. 




(25) 
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4.1.2. Correlations: Having witnessed a good description of charm production as well 
as J/ ip production, we now move to the main topic of the present work. In the following 
we give our results for azimuthal, rapidity-difference, transverse momentum, and jet- 
radius correlation for charm and bottom quarks at 2.76 and 5.5 TeV for pp collisions. 
Deviations from these would signal medium modifications in case of nucleus-nucleus 
collisions. 

Fig. |6] shows pt and rapidity integrated A0 distribution for heavy quarks at a/s 
= 2.76 TeV, 5.5 TeV and 14.0 TeV for both leading order and next-to-leading order 
calculations. As expected the contribution rises with the energy available in the centre- 
of-mass system. It is felt that if our argument about heavy quarks not changing direction 
of their motion due to soft collisions with partons is valid, then drag (or energy loss) 
alone will not drastically alter this feature. It is needless to repeat that at LO all the 
heavy quarks will be produced back-to-back resulting in a peak at A0 = vr. However, 
if the heavy quarks thermalize and flow with the medium, this picture may undergo 
change. We shall come back to this. 

In Fig. [7] we show our results for the transverse momentum, rapidity, and invariant 
mass distribution of charm and bottom quark pairs produced in pp collisions at a/s=2.76 
and 5.5 TeV. Recall that the pair momentum will be balanced by the momentum of the 
recoiling parton. Thus tagging on a high transverse momentum recoiling parton in 
the case of heavy ion collisions can give interesting details of how heavy quarks and 
(mostly) gluons behave in the medium produced in the collision. These results also 
contain a very interesting situation. Consider a heavy-quark produced in LO pQCD in 
a nucleus-nucleus collision. They will be produced back-to-back and are most likely to 
cover different part and length of the system, before they fragment (or coalesce with a 
light quark) to form a D-meson. Thus they would lose a differing amount of energy and 
acquire a net-transverse momentum which was initially identically zero. At least the 
co-linear heavy-quarks produced during splitting of a off-shell gluon would, on the other 
hand, cover similar distances under similar conditions in the plasma, and thus their net 
transverse momentum will remain largely unaltered. It would be interesting to study 
such cases in future more detailed experiments. 

We show our results for rapidity correlation where. Ay = — 7/2, of heavy quarks 
produced in such collisions in Fig. [HI We note that this correlation peaks at vanishing 
rapidity difference. We have also given the LO results for this along with a scaling of the 
LO results with a factor cr^LO/^LO demonstrate that the NLO results can not, in 
general, be approximated by a i^' factor multiplying the LO results, and the inadequacy 
of this shows up most strongly near Ay equal to zero. It is also likely that the rapidity 
difference, specially when the two rapidities have opposite signs may encode effects of 
longitudinal flow in case of nucleus-nucleus collisions. 

Fig. [9] shows the results of our calculation for the jet-radius, i?, correlation, where 
R = A/Arp~+~A0^. It brings out the interesting differences between results for the 
leading order and next-to-leading calculations. Thus, while at leading order we do 
not have any contribution for R < tt, there is a substantial contribution coming from 
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next-to-leading processes for < R < vr. 
4-2. Lead Lead Collisions 

Now we proceed to our results for collision of lead nuclei at 2.76 ATeV and 5.5 ATeV. 
In Fig. [10] we show our results for azimuthal distribution of heavy quarks produced 
from initial (prompt) collision of partons, having transverse momenta of 1-4 GeV and 
rapidities close to zero. The results for LO calculations, having a peak at A0 = vr 
are given, to demonstrate the importance of using NLO results as a base line for these 
studies. We see a sharpening of the coUinear and back-to-back correlations as the 
momenta of the quarks increases, while the correlation, with the exception of the peak 
at A(f) = TT, gets more flat, as NLO processes have a larger role, as the available energy 
increases. We also find less production of pairs of bottom quarks with smaller A0 at 
the same energy, compared to charm quarks, as expected. 

We show our results for production of heavy quarks from multiple scattering of 
jets in Fig. [TTl We have limited our calculations to contributions from quarks and 
gluons having > 2 GeV. Let us first consider our results for charm quarks. At 
both the energies, we see that while charm quarks having transverse momenta around 
1.5 GeV, show a correlation which rises smoothly as we go from collinear to back-to- 
back correlations, the charm-quarks having larger transverse momenta give rise to a flat 
distribution for larger A0. We also note that the contribution of multiple scattering 
of the jets, even though smaller at A(f) ^ tt compared to the contribution of initial 
production, is rather comparable at smaller angular separations. We note that as the 
initial and the final times Tj and Tf appear only as a multiplicative factor ln(r//rj), the 
shape the correlation will remain unaffected by any change in their value. 

The bottom quarks show a very interesting trend. For the lowest momentum 
considered, the bottom-quarks are seen to be produced with a fiat azimuthal correlation, 
while as their momenta increase, the distribution becomes more and more collinear. The 
observation about comparable contributions of multiple scattering of jets and initial 
production at A0 < vr, seen earlier for charm quarks, applies to them as well. 

The results for the angular correlations of heavy quarks produced from the passage 
of jets through QGP are shown in Fig. [121 A very interesting and distinct picture 
emerges for these heavy quarks. We see that these productions are dominated by 
collinear contributions, confirming the nomenclature "jet- conversion" (see Ref. [T2| [21]) 
for them. At small A0 their contribution is similar to that from initial production. The 
corresponding results for bottom quark-pairs show similar trends, but those are an order 
of magnitude smaller than the contribution of initial production. 

And finally the results for the angular correlation of heavy quarks produced from 
scattering of thermalized partons is shown in Fig. [131 Firstly, these contributions are 
smaller by more an order of magnitude than the contributions discussed above. However, 
we still discuss their features as these are quite interesting. The azimuthal correlation 
of charm as well as bottom quark pairs is rather fiat for low transverse momenta but 
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changes steadily to back-to-back at the transverse momentum increases. This, we feel, 
happens as heavy quarks having large transverse momenta can only come from collisions 
of partons having large y/s. This would be possible for partons having almost equal and 
opposite momenta, thus leading to heavy quarks which will be predominantly back-to- 
back. 

In order to get a feeling of the relative contributions of different mechanisms, we 
have shown the correlation of charm quarks for px ^ 2 GeV/c and 4 GeV/c in Fig. [HI 
The first thing we note is that the shapes and relative contributions of the processes 
under consideration remain unchanged with the change in the momentum. Next we note 
that the multiple scattering among the high energy gluons and quarks, termed jet-jet 
interaction, givesi rise to a rather flat azimuthal correlation between the charm quarks, 
which is comparable to the results for prompt production for A0 ^ it. The next large 
contribution is due to the jet-thermal contribution, which is rather fiat for A0e [0 : 2] 
and then falls rapidly. Over this region it is again comparable to the prompt and the 
jet-thermal contributions. The thermal contributions are peaked toward Acp = tt and 
are rather small. 

Recall that we have used a formation time of the plasma as 0.1 fm/c, inspired by 
the parton saturation model. A larger value for will leave the jet-jet contribution 
essentially unchanged, as we discussed earlier. However the jet-thermal and thermal 
contributions are expected to drop if the initial time is increased. Thus recalling 
our results from Ref. [10], we estimate that raising the Tj to 0.5 fm/c the jet-thermal 
contribution may decrease by a factor of 2, while the thermal contribution will come 
down by a factor of about 4. 



4.3. Effect of Flow 

We have suggested earlier that the effect of drag or energy loss of heavy quarks alone 
may not be enough to change their direction of motion, and thus the p^-integrated 
azimuthal correlations discussed in this work may not be affected by the energy loss. It 
may change for a given pt due to migration of quarks to the regions of lower px and the 
Pt dependence of the heavy quark production. The flow of the medium can, however, 
affect the angular correlation considerably, if it is large and if the heavy quarks are 
thermalized. In order to estimate the effect of the flow on the correlation of the heavy 
quarks, we use a toy model used earlier by Cuautle and Paic [38] , and more recently in 
Ref. [39], for studying correlations. 

In order to do this, we proceed as follows. We first give a random orientation to the 
quark-pairs from the NLO pQCD calculations (the NLO MNR code, e.g., at LO gives 
pairs with p^-^ = ^^.2=0). Then we place them at {x,y), randomly chosen according to 
the probability: 

f f dxdyTAix,y,b = 0)TB{x,y,b = 0) 
where Tj is the transverse density profile of the nucleus i assumed to have a uniform 
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density of radius R, and Tpj^{b = 0) is the nuclear thickness for impact parameter, 
6 = 0. Assuming a flow, directed away from the centre, we add the flow momentum 
PTf= PTf (r/r) to the momentum of the heavy quark px- 
We use the blast-model [IQ] to write pTf as 

PTf = iPrmq , (27) 

where 

/3. = /3sX 1^-^)' . (28) 

and r = y/x'^ + y'^. We give results for Ps = 0, 0.3, and 0.6. We show our results 
Fig. [15] for two ranges of pt of the charm quarks, < 4 GeV and > 4 GeV. We see 
that even though the azimuthal correlation is more strongly affected for charm quarks 
having lower transverse momenta for reasonable values of the flow, the basic nature of 
the correlation function remains unchanged. It is likely that if the charm quarks are 
not completely thermalized, the effective flow velocity for them could be smaller, and 
then the above observation becomes even more relevant. Note that large values of /3s 
are normally reached only in the hadronic phase. 

5. Summary 

We have calculated azimuthal, rapidity difference, and transverse momentum 
correlations of heavy quark pairs produced in pp collisions at several energies relevant 
for experiments being done at the Large Hadron Collider, using NLO pQCD. Where- 
ever possible, we have discussed how these could change due to final state effects in 
nucleus-nucleus collisions. These results will act as a base-line for similar studies in 
the case of Pb + Pb collisions at the corresponding centre of mass energies/nucleon, 
to determine medium modifications. We have noted that this picture is enriched (or 
complicated) by multiple collisions among the partons having high energy, which can 
give very different correlations of a magnitude comparable to that of initial productions 
considered above. We have argued, but it remains to be verified, that these correlations 
may not be drastically altered due to the energy loss suffered by heavy quarks, as they 
may not change the direction of their motion substantially, due to soft scatterings. 
These may however, be affected by a strong flow of the medium, if the heavy quarks are 
thermalized. 

In a forth-coming publication we shall address the issue of consequences of energy 
loss on these correlations. 
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